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Optical clocks based on ensembles of trapped ions promise record frequency accuracy with good 
short-term stability. Most suitable ion species lack closed transitions, so the clock signal must be 
read out indirectly by transferring the quantum state of the clock ions to cotrapped logic ions of 
a different species. Existing methods of quantum logic readout require a linear overhead in either 
time or the number of logic ions. Here we describe a quantum algorithmic readout whose overhead 
scales logarithmically with the number of clock ions in both of these respects. The scheme allows a 
quantum nondemolition readout of the number of excited clock ions using a single multispecies gate 
operation which can also be used in other areas of ion trap technology such as quantum information 
processing, quantum simulations, metrology, and precision spectroscopy. 


Tremendous progress has recently been made in opti¬ 
cal frequency metrology [1, 2], Optical clocks now reach 
fractional frequency inaccuracies and instabilities in the 
10~ 18 regime [3-8], outperforming Cs fountain clocks by 
two orders of magnitude and vying to serve as a new 
definition of the SI second [9, 10]. Among the promis¬ 
ing candidates for such a redefinition are ion-based fre¬ 
quency standards featuring very small systematic fre¬ 
quency shifts. However, the poor signal-to-noise ratio of 
single-ion systems entails averaging times of many weeks 
to reach a fractional uncertainty of 10~ 18 [11, 12]. Clocks 
based on strings of ions confined in a linear trap promise 
to overcome this limitation [13, 14], Because of unavoid¬ 
able electric held gradients in such a trap, suitable clock 
ion species may have only negligible electric quadrupole 
moments to avoid systematic frequency shifts [15]. This 
requirement is met by group 13 ion species [16], as well 
as some highly charged ions [17-20]. Most of these candi¬ 
dates lack a suitable transition for laser cooling and state 
detection, so that quantum logic spectroscopy [21] is re¬ 
quired for readout. In quantum logic spectroscopy, the 
internal state of the clock ion is transferred by a series of 
laser pulses onto a logic ion of a different species, whose 
internal state can be efficiently detected. However, ex¬ 
isting methods for quantum logic readout require a large 
overhead in either time or logic ions when an ensemble 
of clock ions is used. 

Here, we suggest using a quantum algorithmic read¬ 
out, cf. Fig. 1, which maps the clock signal (that is, the 
number of clock ions left in the excited state by the spec¬ 
troscopy pulses) on a string of logic ions in binary rep¬ 
resentation and therefore requires a number of logic ions 
scaling logarithmically with the number of clock ions. 
Our approach is based on an algorithm originally intro¬ 
duced in the context of entanglement concentration pro¬ 
tocols [22, 23], which can be implemented efficiently with 
a number of elementary gate operations scaling likewise 
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FIG. 1. (a) In an atomic clock based on Nc trapped ions 
the frequency of a laser is locked to an atomic transition e.g., 
in a Ramsey sequence of two tt/2 pulses (Hadamard gates H ) 
enclosing a free evolution time T. Measurement of the number 
N of clock ions in |1) yields an error signal for the deviation 
of the laser from resonance, (b) For ion species lacking the 
cycling transition needed for direct state detection, a quantum 
algorithmic readout can be used to map N onto AT cotrapped 
logic ions of a different species whose state can be detected 
efficiently. 


logarithmically with the number of clock ions. Moreover, 
we show that the algorithm can be executed by means 
of the multi-ion Mplmer-Sprensen (MS) gate [24-27], a 
well-established tool for quantum control of ion crystals 
[28, 29]. In particular, we demonstrate that only a single 
application of a MS gate involving both logic and clock 
ions is necessary. Such multispecies multi-ion gate opera¬ 
tions are under intense experimental investigation at the 
moment [30], and constitute the most costly resource in 
the current context due to the complexity of the normal 
mode structure in multispecies ion crystals. We discuss 
the feasibility of our approach in a case study of 3 Al + 
and 2 Ca + ions taking into account the full spectrum of 
motional normal modes. 
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The suggested readout method actually realizes a 
quantum nondemolition (QND) measurement on the 
clock ions, and therefore opens up the possibility to de¬ 
vise more sophisticated clock protocols where subsets of 
clock ions are interrogated repeatedly. Beyond frequency 
metrology the QND character of our readout technique 
makes it directly applicable to syndrome measurements 
and error correction schemes for quantum computations 
and simulations with single- or multispecies ion crystals 
[31]. In a more general perspective, this work contributes 
to the development of hybrid systems striving to combine 
the complementary strengths of different physical com¬ 
ponents for quantum information processing, quantum 
simulations and quantum metrology [32, 33]. 

Working principle of ion clocks with direct readout — 
We consider a string of TVc clock ions with a narrow-band 
optical transition of frequency uiq between two internal 
states |0) and |1) which provides the frequency reference 
for the clock. The goal is to stabilize to this transition 
frequency wo a laser field of frequency w with an (un¬ 
known) frequency offset A = to — ujq. To this end, pulses 
of light from the laser drive the clock ions, transferring 
them from |0) to the excited state |1) with a frequency- 
dependent probability. In the simplest schemes, such as 
Rabi and Ramsey interrogation, this probability is inde¬ 
pendent for each ion. The clock readout then consists of 
measuring the number of excited ions TV = |l)i(11 

and using it to infer the excitation probability and hence 
the detuning A, which can then be corrected [34]; cf. 
Fig. 1(a). As explained above, a direct measurement of 
the excited state population TV is impractical for many 
interesting species of clock ion. Instead, one can map N 
onto an ensemble of TVl cotrapped logic ions which can 
be detected efficiently, as illustrated in Fig. 1(b). In di¬ 
rect extension of established readout techniques based on 
quantum logic [21, 35], one could use as many logic ions 
as clock ions (TVl = TVc), performing successive quantum 
gate operations to transfer the state of each clock ion to 
the corresponding logic ion. Alternatively, one could use 
a single logic ion, transferring and measuring the state 
of each clock ion in turn. This imposes a long read¬ 
out time as the ion crystal must be cooled between each 
measurement and subsequent state swap operation. Both 
strategies have prohibitive overhead in additional ions or 
number of gate operations and readout time. The last is 
crucial because time spent on readout adds to the clock 
cycle’s dead time and, thus, through the Dick effect, to 
the clock instability [36, 37]. 

Quantum algorithmic readout — From the perspec¬ 
tive of quantum information theory the quantity of in¬ 
terest - the number of clock ions in state |1) - is the 
Hamming weight of the string of TVc quantum bits (a 
number with Nq + 1 possible values between 0 and TVc). 
In the context of entanglement concentration protocols, 
a quantum algorithm has been developed for the indi¬ 
rect determination of the Hamming weight of a quan¬ 


tum bit string [22, 23]. The algorithm uses an ancil¬ 
lary string of TVl = [log 2 (TVc + 1)] quantum bits on 
which the Hamming weight of the TVc primary quan¬ 
tum bits is stored in binary representation. Thus, the 
necessary number TVl of logic ions (ancillary quantum 
bits) scales logarithmically with the number of clock 
ions. Suitable combinations of clock and logic ion num¬ 
bers (TVc,TVl) are, for example, (3,2), (7,3), and (15,4). 
The algorithmic readout is shown in a circuit diagram 
for the simple case (TVc,TV L ) = (3,2) in Fig. 2(a), and 
for the general case in the appendix. Given TVc quan¬ 
tum bits (clock ions) in any state |TV)c with Hamming 
weight TV and TVl ancillary bits (logic ions) initialized 
in 100 ... 0)l, the algorithm effects a unitary transforma¬ 
tion U such that T7|TV)c|00 • • • 0 )l = |TV)c|*i *2 ■ • ■ *at l )l, 
where the bit string i\i .2 ■ ■ ■i n l (in = 0,1) gives the bi¬ 
nary representation of the Hamming weight of the state 
of clock ions, TV = J2n =l 2 n ~ 1 i n . We will denote by 
\N) l = 1*1*2 *jv l )l the state of logic ions represent¬ 
ing TV. The state of the clock ions after one clock cycle 
is a superposition \ip(A))c = Ejv=o ctv(A)|TV)c, where 
the dependence of the amplitudes cjv(A) on the detuning 
carries the clock signal [38]. Application of the readout 
algorithm generates an entangled state of clock and logic 
ions: 

N c 

l'F(A)) =Z/|V>(A))c|0)l = ^ c n (A)\N) c \N) l . 

N =0 

From a measurement in the {|0), 11)} basis of each logic 
ion one can extract the observable corresponding to the 
estimated Hamming weight of the clock ions, TV est = 
E«ii 2rl_1 l 1 )n(l| = YlN% N \^-)h{KV In a perfect im¬ 
plementation the measurement of TV es t on the logic ions 
exhibits exactly the same statistics as measuring TV on 
the clock ions directly. In particular, the error signal can 
be extracted from TV(A) = (\k(A)|TV est |4'(A)) and used 
to correct lo exactly as in a direct readout. 

Implementation based on M0lmer-S0rensen gates — 
Following Refs. [22, 23], the operation U can be im¬ 
plemented by a sequence of quantum Fourier transforms 
and controlled-phase gates as shown in Fig. 2(a). This 
requires a number of gates linear in TVc, and so provides 
little advantage over the complete state swap mentioned 
above (apart from the reduction in the number of logic 
ions). Fortunately, the algorithm can be decomposed 
much more efficiently using multi-ion Mplmer-Sprensen 
(MS) gates [24-26]. MS gates involve driving the ions 
simultaneously with bichromatic laser fields on red and 
blue sideband transitions to collective modes of vibra¬ 
tions for a duration r, and generate unitary transforma¬ 
tions Ums = exp[-T(S'cc + <Sll + <Scl)] on the string of 
TVc + TV L two-level systems where [39] 

N a Np p a/ 3 

i 1 ) 
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FIG. 2. Quantum algorithmic readout illustrated for Nc = 3 clock and N l = 2 logic ions: after a Ramsey sequency (blue, dash- 
dotted box) the clock ions are in a superposition of states |iV}c, where N £ [0,iVc] is the number of ions in excited states |1) 
(their Hamming weight). A quantum algorithmic readout U (red, dashed box) maps N in binary representation onto a number 
AT = [log 2 (IVc + 1)] of logic ions. Detection of the logic ions in the {|0), 11)} basis provides the number N. (a) Quantum 
algorithm for the indirect measurement of the Hamming weight, taken from Refs. [22, 23]. (1) Logic ions are initialized in 
|0) L = |00)l (states |IV)l denote the binary representation of N). (2) The state is quantum Fourier transformed into |0)l- (3) 
Controlled-phase gates rotate the state in Fourier space to |N)l for a state | A)c of clock ions. (4) An inverse Fourier transform 
yields the state |IV}l of logic ions, (b) The same algorithm decomposed in terms of multi-ion Mplmer-Sprensen (MS) gates. 
The Hadamard gates shown in gray need not be executed. Removing Hadamard gates displayed with dashed boxes merges 
the algorithmic readout with the Ramsey sequence, and requires the laser fields in the MS gates to be phase coherent with the 
Hadamard gates in the Ramsey sequence. Explicit, general forms of all gates are given in the appendix. 


Here, crjA denotes the Pauli x operator for ion ( a,i ), 
where a = C, L and i = 1,..., N a , and the frequencies 
r°f and A“/ are defined through 

r T 1 N l + N c k k 

r t= \ dtn ai (t)n 0j (t), -— = £ -4^- ( 2 ) 

We assume that the Rabi frequencies can be ad¬ 

justed separately for each ion in transverse illumination, 
and may also be tuned in time. In the definition of Afj 
in Eq. (2) the are Lamb-Dicke factors and Sk detun¬ 
ings with respect to the k th transverse normal mode 
of (angular) frequency Vk] cf. Fig. 3(a). The detunings 
have to satisfy Sk and we suppose <5i > 0. A 

decomposition of the desired transformation U in terms 
of multi-ion MS gates and single-ion rotations is shown 
in Fig. 2(b). Remarkably, only a single two-species MS 
gate is required when the interaction between clock and 
logic ions, described by S'cl in Eq. (1), is tuned such that 

rg L = 7r2-b' +2) Ag L ( 3 ) 

where i = 1,..., Nc, j = 1,..., N l [40]. Eq. (3) ensures 
that the MS gate executes the controlled phase gates be¬ 
tween logic and clock ions as shown in Fig. 2(a). Read as 
a matrix equation, it expresses Nq x constraints where 
the right-hand side is fixed for a given configuration of 
the ion string and carrier detuning S = Sk + Vk, while 
the left-hand side has to be adapted accordingly through 
the correct choice of Rabi frequencies and the gate 

duration r. In addition to the interactions between clock 
and logic ions (.S'cl), the gate also induces interactions 
among clock ions (Sec) and logic ions (Sll) only. The in¬ 


teractions on the clock ions have no influence on the read¬ 
out results and can therefore be neglected. On the other 
hand, interactions between logic ions have to be compen¬ 
sated using IVl — 1 gates that involve only logic ions, as 
shown in the appendix. The inverse Fourier transform 
on the string of logic ions involves N l — 1 single-species 
MS gates giving a total number 21 Vl — 1 of multiqubit 
gates (see Appendix) that grows logarithmically with the 
number of clock ions Nq. 

The remaining challenge is to find experimental pa¬ 
rameters that satisfy Eq. (3) . If each ion can be ad¬ 
dressed individually, the time dependence of can 

be exploited to exactly fulfill Eq. (3); as shown in the 
appendix. Here, we will instead focus on finding an 
approximate solution to condition (3) using only time- 
independent Rabi frequencies. This restriction substan¬ 
tially simplifies the experimental implementation, which 
can then employ a single laser beam for each species with 
a tailored intensity profile, rather than requiring full in¬ 
dividual addressing with independent pulse shaping for 
each ion. In this case the left-hand side of Eq. (3) be¬ 
comes r)j L = Qci^Lj'r, which is patently a rank-one ma¬ 
trix while the rhs of Eq. (3) generically has full rank ( N l). 
However, if the detuning <5i is small on the scale of the 
mode spacing Au = v\ — v^, the ion interactions in the 
MS gate are mostly mediated via the highest-frequency 
collective mode. In this regime the rhs of Eq. (3) will 
be almost of rank one in the sense that its singular value 
decomposition exhibits one dominating singular value. 
Accordingly, the rhs of Eq. (3) can be well approximated 
in terms of a rank-one matrix by determining its singu¬ 
lar value decomposition and setting all but the largest 
singular value to zero as detailed in the appendix. We 
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(a) (b) 1 



FIG. 3. (a) Schematic level diagram with relevant frequencies 
during MS gates: atomic transition wo, bichromatic driving 
fields at ujo±8, detuning 5k from sideband Vk due to collective 
mode k. (b) Average readout error probability P elI (red, solid 
line) and multispecies MS gate duration r (blue, dashed line) 
versus detuning <5i. 

emphasize that this is not equivalent to simply dropping 
lower-lying modes [which would correspond to restricting 
the summation in Eq. (2) to k = 1]. 

Case study with Al + and Ca + ions — In the remain¬ 
der of this Letter we discuss in more detail the realization 
of multispecies MS gates with constant Rabi frequencies 
for the concrete case of three Al + clock ions and two 
Ca + logic ions. We assume the ions are held in a linear 
rf Paul trap with soft confinement along the crystal axis 
(z axis) and much tighter confinements in the transverse 
directions that effectively restrict oscillations to only one 
direction (x axis). An important requirement for the 
MS gates is truly collective normal modes; i.e., each ion 
must be involved in the motion to interact with the oth¬ 
ers. For the case of multispecies crystals, these modes are 
found in a stable regime close to the critical confinement 
where the linear chain changes to a zigzag configuration 
[31, 41, 42]. In the appendix we show for the case of 
Al + -Ca + crystals that for moderate trap asymmetry the 
highest-frequency mode is sufficiently collective, and ex¬ 
hibits Lamb-Dicke factors on the order of 10 -1 for clock 
ions and 10 -3 — 10 -2 for logic ions while exhibiting a 
substantial frequency gap Av ~ 27t 500 kHz to the next 
mode. The approximate solutions to Eq. (3) will be ap¬ 
plicable when the detuning in the multispecies MS gate 
is small enough, that is, <5i <C Av. For a given detun¬ 
ing <5i, we require a gate duration r = 2-7r/<5i in order to 
suppress spurious effects of off-resonant terms during the 
MS gate, which implies a time scale in the sub-ms regime 
for the gate. 

In order to assess the performance of our readout strat¬ 
egy, we evaluate the average error probability defined as 
Pe' rr = 1 - EnIoPn P [JV|pc(JV)]. Here P [iV| PC7 (iV)] is 
the conditional probability to correctly receive the Ham¬ 
ming weight N as a result of the quantum algorithmic 
readout given that exactly N ions were in the excited 
state in the initial state of clock ions pc(N). The pn de¬ 


note the a priori probabilities for the states pc(N). In 
the context of ion clocks, the maximum error signal will 
be attained when each clock ion is in an equal superpo¬ 
sition of |0) and 11), which implies that pn = ( I ^)/2 Nc 
and pc{N) is the symmetric Dicke state with N clock 
ions in |1). The conditional probabilities are evaluated as 
P[N\ PC (N)] = tr C L{\N) L (N\Upc(N)® |00) L (00|Wt}, 
where U is the circuit shown in Fig. 2(b). For a given 
detuning the Rabi frequencies in the single multi- 
species MS gate follow from the vectors correspond¬ 
ing to the largest singular value of the rhs of Eq. (3), 
as explained above. Figure 3 shows the overall er¬ 
ror probability P err and the gate duration r versus de¬ 
tuning <5i in the MS gate with otherwise ideal pro¬ 
cesses. For example, a detuning <!>! = 2-7r20kHz cor¬ 
responds to a gate duration of r = 50 ps and an er¬ 
ror probability of 0.5%, which is achieved with Rabi 
frequencies (Qli, ^ 2 ) = 27r (730.34, 363.81) kHz and 
(H C i, Hc 2 , H C3 ) = 2tt (49.71, 45.84, 53.81) kHz for logic 
and clock ions, respectively. In this example the carrier 
transition is driven off resonantly with a Rabi frequency 
to detuning ratio of at most H a ,/<5 = 0.23. However, 
the associated ac Stark shifts cancel due to the bichro¬ 
matic, red- and blue-detuned drive and can therefore be 
neglected. At the same time, the 50 ps gate duration is 
short enough compared to the 1.17 s D 5 / 2 state lifetime in 
Ca + that readout errors due to spontaneous emission will 
be on the per mill level. The Rabi frequencies required in 
the present example can be tailored by implementing the 
gate with a tightly focused TEMi 0 laser beam addressing 
the three clock ions, such that each ion is located at a 
position in the transverse intensity profile corresponding 
to the Rabi frequency given above. 

The more complex spatial structure of Rabi frequencies 
needed in a longer string of ions can be engineered with 
spatial light modulators [43] or multichannel acousto- 
optical modulators. Combined with the freedom to per¬ 
mute clock and logic ions and to choose different solutions 
to Eq. (3), we expect experimentally feasible implemen¬ 
tations for more than 15 clock and 4 logic ions. Impor¬ 
tantly, readout errors will affect only the stability, not 
the accuracy, of the ion clock. 

Conclusion — The quantum algorithmic readout sug¬ 
gested here performs a QND measurement of the Ham¬ 
ming weight of clock ions and may therefore allow 
more complex clock protocols using repeated readouts of 
(sub)ensembles of ions or preparation of the clock ions in 
Dicke states for nonclassical frequency metrology. We en¬ 
vision that the techniques developed here in the context 
of atomic clocks may prove useful also in other areas of 
ion trap technology, such as quantum information pro¬ 
cessing, quantum simulations, metrology, and precision 
spectroscopy. 
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Appendix 


A. Calculation of normal mode spectrum for transversal oscillations in a multi-species ion crystals 

For sake of completeness we summarize the determination of the normal mode spectrum for the transversal os¬ 
cillations of a two-species ion chain, following the treartment of [42]. In the example setup described in the main 
text two logic ions (Ca + , mass = 40amu, laser wavelength A^ = 729.1 nm) and three clock ions (Al + , mass 
me = 27 amu, laser wavelength Ac = 267.4 nm) are trapped along the crystal axis (z-axis). The trap frequency 
for the logic ions is = 27r874kHz. The asymmetry parameters are defined as a = v x /vz and o. yx = v J y /v k . 
We fix a yx = 5 to suppress oscillations along the y-axis. The normal modes are calculated for different values of 
a, because the oscillations in ^-directions will be used for the gate. As sketched in Fig. 4(b) the clock ions (index 
k = 2,3,4) are placed in the middle and the logic ions (index k = 1 and k = 5) on the outside. With masses 
(mi, m 2 , m 3 , 7714,7715) = (mi,mc,mc,mc,mi) the kinetic energy is given by 

5 

T (p) = J2p*/ 2rnk - 

k =1 

The potential energy due to the electrostatic component of the trap is the same for each ion, as all ions carry a single 
positive elementary charge q. Assuming a radially symmetric electrostatic trap the total electrostatic potential is [41] 

5 

Vs(£) = hYl ( b ° z k ~ I b ° x k ~ vl) > 

fc=1 

The potential strenth is parametrized by the parameter bo, in units of energy divided by length squared. To keep 
the ions also radially trapped, an additional time-dependent radiofrequency potential is used to create an effective 
mass-dependent potential in x and y direction 


V RF {x) 


IJ2 ~ ( bxXk+byV '*) • 

k =i mk 


The potential strength is parametrized by b x and b y , in the same units as bo. The parameters bo, b x and b y are fully 
determined by the physical parameters of the setup: logic ion mass mj,, logic ion frequency vjf and the asymmetry 
parameters a and a xy . It is important to note that the two species of ions experience different radial potentials because 
the pseudopotential generated by the radial AC fields in a Paul trap is mass-dependent. As a result, lighter ions feel 
a tighter transverse potential. The corresponding trap frequencies of clock and logic ions are denoted by u x and iA, 
where v\: < v x since WAi/mca = 27/40. In addition to the trapping potential the ions interact via the Coulomb 
repulsion potential 


Vi(x) = £ 

k>j 


47T£o 



With these definitions we can write the total energy of the system as E(x,p) = T(p) + V{x), where V(x) = Vs(x) + 
V R f(x) + Vi(x). 

We find the steady state position x 0 of the ions by numerically minimizing [45] the potential energy V under the 
condition z\ < < z^. In this study we choose large enough asymmetry parameters a so that the solution is 

always a linear chain without zigzag configuration [31, 41, 42]. As the oscillations around the steady state will be 
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small, we use second order Taylor expansion to obtain an approximate harmonic potential. The different directions 


x, y and 2 decouple in this approximation. Denoting p *. = m±±k and Vjy = d XkX V(x) 
^-direction only is 

5 


the energy for motion in 


P 2 1 

E x — ^2 + X! nVkj x k x j- 


1, 


fc=i 


2 m fc 


k,j 


In coordinates with scaled position Xk = yj and momentum pk = normalized to mass mo = 1 amu, the 

kinetic term becomes diagonal and the potential transforms as V^j = y mk m Vkj■ In these coordinates E x reads 


p2 L _ 

E* = y: + 5Z b ^kjXkXj- 


fc=1 


2mo 


k,j 


For the normal modes we numerically diagonalize V = ODO T with a dimensionless orthogonal matrix O and diagonal 
matrix D of dimension frequency squared times mass. The eigenfrequencies are then given by 

»k = yjDkk/mo. 


In analogy to [42] the Lamb-Dicke factors of a mode k for an individual ion with index j are 


k 27T 

n = —O 
h ^o 



with A j the wavelength of the laser addressing the j -th ion and Ojk the j-th entry of the eigenvector for the k -th 
normal mode. 


(a) normal mode frequencies 



(b) 

KT 1 

77 ai 

nr 2 
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Lamb-Dicke factors 
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0 ^ 

LI Cl 

0 0 (0 

C2 C3 L2 




__PL1 = PL2 


FIG. 4. (a) Normal mode frequencies of transverse vibration for a crystal of Nc = 3 Al + clock ions between jVj, = 2 Ca + 
logic ions (see b, inset) versus asymmetry parameter a = (ratio of transverse to axial trap frequencies for logic ions, 

= 27t 874 kHz). (b) LD factors and 77 jh of, respectively, the three clock and two logic ions for the normal mode of highest 
frequency (thick line in (a)). For large asymmetry parameter transverse motions of clock and logic ions decouple due to the 
mass-dependent transverse confinement. The dotted line in (a) and (b) marks the case studied in the text. 


In Fig. 4(a) we show the spectrum of the transverse normal modes for the case of (Nc, IVl) = (3, 2) in the ordering 
shown in the inset of Fig. 4(b). The generic result is that for large asymmetry parameter a the Ay + Nq normal modes 
split into two groups involving either the AL logic or the Nq clock ions. To have truly collective normal modes involving 
both species of ions (as required for the Mplmer-Sprensen gate) the asymmetry must be kept moderate. However, 
below a critical value of a the normal modes become unstable and the crystal changes to a zig-zag configuration 
[31, 41, 42]. For moderate trap asymmetry the highest lying modes is sufficiently collective, and exhibits Lamb-Dicke 
factors on the order of 10 -1 for clock ions and 10 -3 to 10~ 2 for logic ions while exhibiting a substantial frequency gap 
Au to the next mode. For an axial trap frequency of 27r874kHz and an asymmetry parameter a = 2.5 the highest- 
lying transverse mode has a frequency = 2tt 3.14 MHz and the gap to the next mode is Av = zn — = 27r480kHz. 






















LI 

Cl 

C2 

C3 

L2 

Vi 

= 2tt3.14 MHz 

0.007 

0.098 

0.113 

0.098 

0.007 

V2 

= 27t 2.66 MHz 

-1.301 ■ 10 -2 

-0.133 

-2.429 • 10" 8 

0.133 

1.301 - 10 -2 


= 2tt2.06 MHz 

-0.024 

-0.060 

0.137 

-0.060 

-0.024 

Vl 

= 2tt 1.79 MHz 

3.996- 10~ 2 

-4.318 • 10 -2 

-1.079- 10" 6 4.318- 10~ 2 

-3.996 ■ 10" 2 

V5 

= 2t r 1.72 MHz 

0.034 

-0.067 

0.071 

-0.067 

0.034 


TABLE I. Numerical solutions of normal mode frequencies and LD factors r/ ai for all ions and modes at a = 2.5. Each row 
corresponds to a normal mode as indicated by the index of the corresponding frequency. The symmetry in LD factors is due 
to the symmetry of the ion configuration. 


Assuming all laser beams are aligned with the rr-axis, the corresponding LD factors for logic and clock ions are given 
in Tab. I. The LD factors for the highest frequency mode are most relevant for the gate dynamics. Slightly larger LD 
factors were used in [46, 47] to achieve high-fidelity Mplmer-Sprensen (MS) gates on two ions and entanglement of 
up to four ions, respectively. Other demonstrations of high-fidelity multi-ion gates used LD parameters which were a 
factor of two to three smaller than the values considered here [48-50] . The information on the normal modes and the 
gate detuning is enough to calculate the r.h.s of Eq. (3) and find appropriate Rabi frequencies. 


B. Derivation of the effective Mplmer-Sprensen-Hamiltonian 

A Mplmer-Sprensen (MS) gate is achieved by the interaction of multiple ions with two laser fields equally detuned 
to the upper and lower sideband of a collective motional mode [24, 27]. The two lasers A and B (red and blue detuned, 
respectively) have frequencies u ja = loq — 6 and w B = wo + d where 5 is the detuning to the carrier frequency ujq. 
Therefore 5k = 5 — gives the detuning with respect to the motional mode v k mediating the interactions. Here <5i > 0 
is supposed to be a small detuning to the motional sideband with highest frequency in Fig. 4(a). First we assume 
a general setup of N arbitrary ions with individual Lamb-Dicke factors t]a/b j f° r i° n 7 and collective mode k. The 
bichromatic laser has, in general, time dependent Rabi frequencies Q A/B,j(t ) and a laser phase fa that is assumed to 
be equal for both lasers A and B but possibly different for each ion. The Hamiltonian of the system consists of the 
internal energies, motional energy and the interaction with the lasers. Changing into an interaction picture gives a 
time dependent Hamiltonian in Lamb-Dicke expansion (h = 1) 

N 

U(t) = ^2 nA,j(t ) t/aj (e^e^'crt a k + h.c.) + n Bj -(t) (e~ lSkt e~ l<t>j cr+a f k + h.c.) . (4) 

j,k=l 

Here a\ and a k are the creation and annihilation operators of the collective motional mode with frequency The 
Lamb-Dicke approximation was used keeping terms only to linear order in il k A / B j an d a rotating wave approximation 
was applied, neglecting all terms rotating faster than 5jy. The laser phases can be absorbed through a unitary 
transformation 77(f) = VHV ^ where 



such that 


N 

U(t) = ^2 toA,j(t ) 11a, j (e lSkt afa k + h.c.) + 0 BJ (i) (f 




aja' k + h.c.) . 


j,k=l 


The effective Hamiltonian can be calculated from TL{t) using a Magnus expansion 

U (At) = exp (Y (At)) 

with 

H ) 2 f At 


r At 


Y(At) = —i ^(t) dt + 


/ o U o 


H(t")dt" ,U(t) 


dt 
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Oscillating integrands proportial to e l ^ s+v ' l ~ Vk ' >t or even faster phases are assumed to average to zero over the gate 
duration and therefore will be neglected. The unitary time evolution is then given by a time dependent effective 
Hamiltonian 

U(At) = exp i H e ff (A t)j . 
and 'H e s can be written in a compact form as 

N 1 f At 

T-L c g (At) « Y, —— / S^. k (; t ) + Sy k (t) + S z>k ( t) dt 

k =1 4 k 

where 

N 

S x ,k ( t ) = Y ( n A,j ( t) r)A,j + &B,j (t) Vbj) fa, 

3 = 1 
N 

Sy,k ( t ) = Y ( t) riAJ - n B,j (t) VB,j) fa 

3—1 
N 

S z ,k (t) = Y 2 (*) Vaj + Q B,j (t) VBj) {^A,j (t) Vaj - ( t ) f/Bj) (2aJ.Ofc + l) ctJ 

j =i 

This representation emphasizes the different contributions to the effective Hamiltonian. 5^ fc and fc give rise to 
the usual collective spin flips in a MS gate and S z>k are energy shifts for the internal states of the ions. Note that 
both S y> k and S Z:k are proportional to the differences in Rabi frequencies and Lamb-Dicke factors for the lasers A 
or B and therefore vanish if those are equal. This Hamiltonian is an approximat solution assuming slowly varying 
Rabi-Frequencies compared to the detuning on the timescale of the gate duration. Now the unitary transformation, 
with V given above, is applied to Ties to find the full Hamiltonian of the MS interaction, namely 

N r*/\t '2 2 

Hus (At) = vn cS (At) yt = Yjt (Xh) (*) + (Xa)" (*) + s z<k (t) dt 

k=1 ^ k J 0 l J 

with the operators 

N 

s£,k (*) = E ( Qa ’3 (*) lAJ + n Bd (t) n B,j) (fa cos fa - fa sin fa) 

3=1 

N 

Sy,k ( t ) = Y (*) dAj ~ n B,j (t) Vbj) (fa cos fa + a* sin fa) 

3=1 

and S z k (t ) stays unchanged. 

If we assume now that each ion interacts with both laser beams in the same way, meaning that f Iaj (t) = H B,j (t) = 
flj (t) and r\ k A • = r)% • = fa hold for each ion j and mode k, the effective Hamiltonian reduces to interactions given 
by the S% k (t) terms only: 


%ms (At) = "Y ^ ^E % (*) fa (fa cos fa ~ fa sin fa) j ^E W fa’ [fa' cos fa> - fa, sin j dt 

Adjusting the Laser phases fa allows to alter between a x a x or a v a y interactions and also determines the overall sign 
of these. Choosing fa = fai = i r simplifies the Hamiltonian to 


N N 1- fc 


Tj .71 . t’AI 

T~Lms (At) = E EhE^ / Ctj (t) Clj’ (t) dt 

3,3' = 1 k= 1 k J ° 
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Finally, to compare this result to the quantum gate used in the main text, we define Tjj/ (At) and —— as 


Aii' 


/•At i n 

Fjj' ( At ) = / % (t) (*) dt and -— = 

J o A 77 fe=l 


N h h 

TO 


This gives 


where 


17ms (At) = exp (-iS (At)) 


AT 


5(At)= £ ^ 

i,i'=i " 


Fjj' (At) an an 

A ., 


This is identical to the unitary evolution used in the readout strategy if we label the N ions accordingly as clock or 
logic-ions, i.e j —> ai. 


C. Algorithmic Readout for General Number of Ions 


Here we explain the readout algorithm - which was discussed for the simple case Nq = 3 , = 2 in the main text 

- for arbitrary number of clock and logic ions Nq and TVl. The generalisation of the schematic circuit in Fig. 2(a) to 
the case of arbitrary numbers of ions (following [22, 23]) is shown in Fig. 5(a). The Ramsey-sequence is performed on 
all clock ions (leaving them in a superposition of states |7V)c), and the logic ions are prepared in the ground state at 
the beginning of the algorithm IA. After applying the quantum fourier transformation on the logic ions the algorithm 
consists of controlled phase gates using the clock ions as control bits. This way an excited clock ion gives an additional 
phase to the excited state in each logic ion via the unitary phase gate 


Rk = 


0 

e 27Ti/2 fc 


and generates the state |IV)c|N)l- The inverse quantum fourier transformation produces the Hamming-weight N in 
binary representation as the logic ions’ state, which can then be retrieved by a measurement on the logic ions. 

Generalising the implementation of the readout by means of Mplmer-Sprensen (MS) gates requires some more work 
than the schematic description, mostly due to the inverse quantum fourier transformation. The initial quantum 
fourier transformation in IA is again given by Hadamard gates on the logic ions, since they were initially prepared 
in the ground state. The main part of the algorithm, i.e. the controlled phase gates between clock- and logic ions, 
are performed by a single MS gate and additional single qubit phase gates. This MS gate is described in the main 
text and Eq. (3) shows the condition that the coefficients need to fulfill in order to implement the desired algorithm. 
Additional to the desired dynamics a general setting will also involve MS interactions between all pairs of logic ions 
with interaction strengths T^/AT;, where j,f = 1 , • • • , IVl- Applying a correction dynamic U COII can reverse these 
interactions. Wherein f7 corr consists of Nj j — 1 Mplmer-Sprensen gates each involving a smaller subset of logic ions 
and the signs of the interactions are inverted by changing the laser phase for a different ion in each gate. In Fig. 5(b) 
the single phase gates associated with 17 ms are labeled as s*. Those are controlled phase gates 


$k — 



(5) 


with phases 9k = Nc 7r 2~ k for k = 1,..., W. Additional single ion phase gates (and Hadamard gates) on the clock 
ions are left out in this discussion as well as in the given figures (displayed in grey). These gates are only necessary 
to recover the initial state of the clock ions, but this is not the focus of the algorithmic readout presented here. 
Hadamard transformations surrounding the MS gates are used to relate the physical a x ® cr x to a z ® a z . 

With more logic ions, the implementation of the inverse quantum fourier transformation become more costly in 
terms of the MS gates needed. An efficient circuit for the quantum fourier transformation consists of blocks with 
different controlled phase gates separated by Hadamard gates, see e.g. [51]. In analogy to other parts of the readout, 
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Ramsey-Sequence Readout LA 



(b) 


(i) 

|w)cio) L 


( 2 ) 

|w)c|0)l 


(3) (4) 

|W)c|N) l |JV> c |JV)l 



|W)o|0)l |W)c|0)l |JV) c |N) l \N)c\Kh 


FIG. 5. Full algorithmic Readout for a general number of clock- and logic ions Nc,Nl. (a) As described in the main text, 
the clock cycle starts with a Ramsey sequence (blue, dash-dotted box) on all Nc clock ions, uses the quantum algorithm U 
(red, dashed box) to encode the number of excited clock ions onto the AT logic ions and gives Hamming weight in binomial 
representation with a subsequent measurement of each logic ion. So the unitary transformation U starts at (1) with the state 
|7V)c|0) l where N clock ions are in the excited state and all logic ions are prepared in the ground state. Then the quantum 
fourier transform is applied to the logic ions to give |./V}c|Q}l as the resulting state at (2). The controlled phase gates Rk add 
additional phases to the logic ions for each excited clock ion to give the state |A r )c|N)L at (3). The inverse quantum fourier 
transform yields the state |-/V)c|A[}l at (4) so that the detection of each logic ion gives the binary representation of N. (b) 
Decomposition into quantum gates for an experimental realisation. In the first step the quantum fourier transformation is 
performed via Hadamard gates on the logic ions. Then a Mplmer-Sprensen-gate connecting clock- and logic ions transfers the 
information about the Hamming weight onto the logic ions and the gates in C/ corr correct for unwanted dynamics. The inverse 
quantum fourier transformation is decomposed into AT — 1 Mplmer-Sprensen-gates and related single ion phase gates separated 
by a Hadamard gate. For simplicity each phase gate is labeled s although they all describe different phase shifts. The detailed 
phases for s and the coefficients for each Ums operation are given in the text. Again, Hadamard gates shown in grey drop in 
pairs and do not need to be executed. Grey dashed gates merge the readout with the Ramsey sequence. 


these steps are performed using MS gates and single ion phase gates. The inverse quantum fourier transformation 
on the logic ions is implemented by a series of Ay — 1 such steps, each involving an increasing number of ions and 
ending with a Hadamard gate on the last ion involved. Every MS gate is described by the same mechanism as given 
by Eq. (1) in the main text. It is therefore determined by F'yb/A^; on the logic ions. For the step involving ions 1 
to k the coefficients are 
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■pLL 

-fp = —7r • 2 J+J '- 3 and 

A LI ; 

n 


pLL 

L jk 


= —7r • 2 J 


-/c-l 


pLL 

with j,j' = 1,2, ...,fc — 1. For k = AA the largest coefficient ~ll 1,JVl 2 = — n • 2 27Vl_6 gives the extreme case, 

AjV L -l,jV L -2 

requiring the largest Rabi frequencies. This can set a limit to possible implementations of this algorithm for large AA- 
The corresponding single ion phase gates in this step are also described by Eq.(5) with phases 


Oj = -27r • 2" (fe - j) 


for j = 1, 2,..., k — 1 and 


k-1 

0 k = -2t r Y 2- (fe " m) 

m =1 

for the &-th ion. These coefficients are chosen such that they give the desired controlled phase gates and also discard 
undesired interactions among the logic ions. 


D. Exact and approximated gate dynamics 


Exact solutions — Exploiting the possible time dependence of £ l a i(t ) leads to exact solutions of the desired dy¬ 
namics. The Rabi frequencies required for this can be obtained by considering the singular value decomposition for 
the r.h.s. of Eq. (3). Let s n > 0 be the corresponding singular values (of dimension Hz) and v n £ IR^ 0 and w n £ IT V| ' 
the left and right singular vectors, respectively. Eq. (3) then is equivalent to 


JVl 


d tn C i(t)^Lj{t) = Y S n V i W j- 


( 6 ) 


Let P n (t) be any set of orthonormal functions on [0, r], i.e. 

[ dtP n (t)P m (t) = 6 nm - 

Jo 

When the Rabi frequencies are chosen as 

W 1 

Pici(t) = Y S n — Vi Pnit) 
n—1 r " 

j{t) =Y TmW 7 P mij ) 

m= 1 

with r n £ R Eq. (6) is fulhlled. Possible choices for P„(t) include piecewise constant functions (in the spirit of 
segmented MS gates explored in [52]) or e.g. Legendre polynomials. A change in the sign of Rabi frequencies may be 
necessary therein and can be reached by a change of 7r in the laser phase on the respective ion. The additional vector 
r can be used to further adapt the required Rabi frequencies. 

Approximation scheme — Alternatively, and much more appropriate for the context of optical ion clocks, an 
approximate solution to condition (3) can be achieved with time-independent Rabi frequencies: Here the left hand 
side of (3) becomes T^ L = Plci^hj T which is a matrix of rank one. The r.h.s. of Eq. (3), which generically has full 
rank (AA), can be approximated by a rank one matrix as follows: If the detuning Ai <C Av = V\ — is small as 
compared to the mode-spacing the MS interaction is mostly mediated via the highest lying collective mode with small 
corrections from other modes. 

When the contribution of lower lying modes is neglected the interaction factorises as 1/A“^ « r] a ir]pj/5i and 
also becomes a rank one matrix (single mode approximation). Therefore the Rabi frequencies fla = <$i/4 ij^ and 
= 2 _ W +1 1 A /r]} :) with i = 1 , • ■ • , Nq and j = 1 , • ■ • , AA give an approximate solution. 
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A better result can be achieved by finding the singular value decompostion of the r.h.s. of Eq. (3), and setting all 
but the largest singular values to zero (rank one approxiamtion ), such that Eq. (3) becomes 


= Si v\w]. 


(7) 


This is solved by setting 


i = 



— 



rw] 


for any convenient choice of r 6 R. In the case study discussed in the main text r was chosen such that the largest 
Rabi frequency Uli still has the same values as in the single mode approximation, that is = ^i/4%i- This ensures 
that the adiabaticity condition in the MS gates is fulfilled for all ions. 


